function phi = eq_phi(x,y,m,n)
phi = e^(1i*2*m*x+1i*2*n*y)/pi;
